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Abstract 

Some interesting optical instruments such as the Eaton lens and the Invisible 
Sphere require singularities of the refractive index for their implementation. We show 
how to transmute those singularities into harmless topological defects in anisotropic 
media without the need for anomalous material properties. 
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1 Introduction 



Recent developments in the field of transformation optics flflo] and metamaterials 10 15 



have led to the proposal of a number of remarkable optical instruments, including the Eaton 
Lens 16 19 , the Invisible Sphere [9,20 21 and invisibility cloaks |H[5| [22| - [24] . However, the 
implementation of these devices often requires physical parameters that are very difficult 
to achieve in practice. For example, a requirement for some is that the speed of light goes 
to zero (e.g. at the geometric centre of the Invisible Sphere and the Eaton Lens) or to 
infinity (e.g. in conformal invisibility cloaking |4||22]). Such extremal requirements can be 
usually met only for discrete resonant frequencies, which severely restricts the applicability 
of these devices [24j . 

A method has been proposed by Tyc and Leonhardt to transmute the refractive index pro- 
files of singular optical devices into equivalent non-singular profiles with a slight anisotropy 
25 . This technique is based on the expansion of space around a singularity with a non- 



constant factor to bring all index values within a finite range. The transmutation of singu- 
larities was a crucial step towards the practical realisation of singular optical devices (e.g. 
the Eaton lens [19]). However, the method of transmutation established in [25] always 
leads to anomalous eigenvalues for the permittivity and permeability tensors (i.e. values 
that are less than unity), which means that the transmuted profiles can be implemented 
only for discrete resonant frequencies 24 . In order to avoid anomalous transmuted index 
values a modification of the original method was suggested by Banner 26 . The modified 



technique is based on changing the expansion factor derived in 25 to make all transmuted 
values greater than unity. 

In this paper we show that a completely non-anomalous transmutation can also be achieved 
by making the transmutation derived in |25j| partial. This modification makes transmu- 
tation a local operation, which can be advantageous in a number of cases. First, partial 
transmutation minimises the anisotropic region in transmuted singular devices. Second, 
this technique can be applied to devices with potentially any number of isolated divergent 
singularities. Third, partial transmutation leaves the major part of the refractive index 
profile intact, which allows us to exploit desirable features of the non-transmuted isotropic 
refractive index profile that would otherwise disappear in a full transmutation done ac- 
cording to (25] or 



26 



For our purposes the third point is the crucial one, which enables 
us to produce a partially transmuted Invisible Sphere, which can be used as an invisible 
background profile for invisibility cloaks to avoid superluminal propagation while achieving 
electromagnetic cloaking. For a full account on the application of the partially transmuted 



Invisible Sphere as a background for invisibility cloaks see 27 



2 



2 Examples 



First, we introduce the idea behind partial transmutation through the examples of the 
Invisible Sphere and the Eaton Lens. The Invisible Sphere is an optical lens that is invisible 
in the limit of geometrical optics [9 ,20 ,2l]. Each incident ray completes a loop in the sphere 
and emerges in the direction of its entry (Fig. lA), creating the illusion that space is empty. 
This device has a spherically symmetric refractive index profile given by 20 



n{r) = [Q - 1/{3Q) ] if r ^ 1, and n(r) = 1 ifr>l, (1) 

where 




Qir) = \h-+\ -. + 7^- (2) 



In the expression above we assume that the device has unit radius. It is easy to see that 
n(r) tends to infinity as r — )■ and that n(r) diverges as ~r~^/^ near the geometric centre 
of the sphere [o] . The diverging refractive index profile implies that the speed of light goes 
to zero at the geometric centre, which is extremely difficult to achieve in practice. However, 



we could follow the method described in 25 to transmute the refractive index profile of the 
Invisible Sphere into an equivalent anisotropic profile by expanding the space around the 
singularity with a factor that increases continuously as we get closer to the singular point. 
Following the transmutation, all eigenvalues of the permittivity and permeability tensors 
will be finite and thus no singularity will be present. Fig. IB shows the light trajectories 
in the transmuted profile. 

However, calculations show that the complete transmutation of the profile of the device 
leads to anomalous eigenvalues for the permittivity and permeability tensors, i.e. to eigen- 
values that are less than unity. This can be a difficulty since such a profile would work only 
for incident waves with frequencies that are close to the resonant frequencies of the imple- 



menting material 28 . In order to avoid this difficulty, we perform a partial transmutation 
on the original profile of the Invisible Sphere. This means that we leave the outer region 
of the profile unaltered and expand space only in the immediate vicinity of the singularity 
within a given radius b. Fig. IC shows the partially transmuted profile. Calculations show 
that if we start the transmutation of the inner region at a radius less than or equal to 
r = 0.2489, all eigenvalues of the permittivity and permeability tensors will be greater 
than unity, yet finite. 

The partially transmuted Invisible Sphere can be used as an invisible background profile 



for anomalous invisibility cloaks 27 , because the cloak can be placed outside the trans- 
muted region, where the refractive index is high. Thus the Invisible Sphere can raise the 
anomalous index values of the cloak above one, thereby making the implementation of the 
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Figure 1: A: Light trajectories in the Invisible Sphere. B: Trajectories in the completely 
transmuted Invisible Sphere. C: Trajectories in the partially transmuted Invisible Sphere, 
where the transmutation starts from r=0.2489. 



cloak practicable over a broad range of frequencies 27 . 

Another important example of singular optical devices is the Eaton lens [7,16,17,19 , which 
is an omnidirectional retroreflector (Fig. 2A). It also has a spherically symmetric refractive 
index distribution given by 



n(r] 



^2/r — 1 if r ^ 1, and n{r) = 1 if r > 1. 



(3) 



In the expression above we assume that the device has unit radius. The refractive index 
goes to infinity as r — >■ 0. It can also be shown that n(r) diverges as ~r~^/^ near the 
geometric centre of the the lens jo]. We can follow the same procedure for the Eaton 
Lens as for the Invisible Sphere to obtain the completely transmuted version of the lens 
which has finite but anomalous refractive index values (Fig. 2B). By performing partial 
transmutation at a radius less than or equal to r = 0.3585, we can arrive to an anisotropic 
profile whose tensor eigenvalues are greater than unity and finite at all points (Fig. 2C). 
In fact, it turns out to be generally true that for a profile with an isolated singularity 
and with asyptotic form n(r) (where — 1 < p < 0) an equivalent partially transmuted 
anisotropic refractive index distribution can be found, where all values of the permittivity 
and permeability tensors are greater than unity and finite, and whose overall influence on 
the light trajectories is equivalent to that of the non-transmuted profile. 



3 Method 

In the previous section we argued that refractive index profiles with isolated divergent 
singularities can be changed into equivalent non-anomalous finite profiles by partial trans- 
mutation. In this section we establish this result in a more general setting. Let us consider 
profiles that have singularities that behave like with —1 < p < in the asymptotic 
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Figure 2: A: Light trajectories in the Eaton Lens. B: Trajectories in the completely trans- 
muted Eaton Lens. C: Trajectories in the partially transmuted Eaton Lens, where trans- 
mutation starts from r=0.3585. 



limit. This implies that the refractive index goes to infinity and the speed of light goes to 



zero as r — )■ 0. It was shown in 25 that we can give enough space for the propagation 
of light and bring all refractive index values into a finite range if we transform the radial 
coordinate according to R = r^"*"^ and leave the angular directions unchanged. The expan- 
sion factor between the virtual and physical space than becomes rjn = dR/dr = {p + l)r^, 
which implies that light will move t^r times faster in physical space than in virtual space. 
Therefore, the refractive index will be = l/(p+l) in the radial direction and ng = = 1 
in the angular directions in physical space. Consequently, the singularity is effectively re- 
moved. 

However, following the transmutation another problem emerges. Since = (dete)^"^ 
(where n and e are the refractive index and permittivity tensors respectively) |9], it follows 
that e = {det n)n~^. Therefore, Er = {p + 1) and = eg = {p + Consequently, for 

— 1 < j9 < the radial component of the permittivity tensor of the profile will always be less 



than one, giving rise to anomalous dispersion 28 , which severely restricts the applicability 
of the device. 

In order to avoid such anomalous refractive indices, instead of transmuting the entire profile 
of the singular optical instrument, we leave the outer region of the profile intact and trans- 
mute only the immediate vicinity of the singularity within a given radius b (Fig. 2C). In 
this case, the expansion of the radial coordinate in the transmuted region is still according 
to the power law: R ~rP+\ but now we will need to include a normalisation factor in front 
R = jf-r^~^^ to make sure that Rib) = b (i.e. that the transformation does not introduce 
impossible jumps in the trajectories), where b is the outer boundary of the transmutation. 
Consequently, the expansion factor will now be rjR = dR/dr = M'^ip + l)rP, giving rise 
to refractive index values = M /{p + 1) and ng = = M . It is easy to show that 
A/'(&) = IP, which tends to infinity as 6 — (Fig. 3). Since n^, ng, and are all propor- 
tional to Af{b), it follows that their value can be increased to any arbitrary finite value by 
choosing a sufficiently small b for the radius of the transmutation and thus the area of the 
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transmutation can be reduced to an arbitrary small area. 



4 Conclusion 

In conclusion, we have established that for optical instruments that have isolated singular- 
ities diverging as (where — 1 < p < 0) an equivalent partially transmuted anisotropic 
refractive index can be found whose tensor eigenvalues are always greater than unity and 
finite, and whose overall influence on the light trajectories is completely equivalent to that 
of the non-transmuted profile. This method can be applied to optical devices with any 
number of singularities. 
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A Appendix 

A.l Material properties 

We calculate the material properties of the transmuted optical devices in physical space 
by the method described in [6||8|[9]. We consider an optical device with a refractive index 
distribution that is spherically symmetric with n(r) = for r ^ 1 and n{r) = 1 for r > 1. 
We transform the radial coordinate by the following function 

Rir) = }, . r^^^ if r ^ 6, and Rir) = r ii r > b. (4) 

We leave the angular coordinates unchanged. This represents a transformation from vir- 
tual space [described by the metric tensor gi'j/ = diag(l, r^, r^sin^^^)] to physical space 
(described by •yjk = diag(l, i?^, i?^sin^^)). The transformation matrix between the two 
spaces IS is given by the Jacobian: A\ = d{r,e,(j))/d{R,d,(f)) = diag(dr/di?, 1, 1). The 
mixed components of the permittivity and permeability tensors fi\ can be calculated 
in physical space by the recipe given by equation (9) in |6j and equation (32.1) in [9] 

^\ = f^\ = n " 9i'j'^'[^''j- (5) 

In the expression above g^^ is the inverse of Qij, g is the determinant of gij, 7 is the de- 
terminant of and n is the refractive index. In the above equation we assume that the 
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dielectric medium is impedance-matched, i-e.e\ = Substituting the relevant expres- 
sions into (6.5) gives 



= = . diag(--,-,-j. (6) 

Substituting n = and R = j^^^ r^"*"^ into the above expression gives the exphcit expres- 
sion for the mixed tensors in physical space: 

e\ = diag(e^, e^, se) = fi\ = diag(AA(6) (p+l), ^) , (7) 

This result shows that the singularity has been effectively removed through the transmu- 
tation and that the tensor components are directly proportional to A/'(6). Therefore, by 
choosing a sufficiently small transmutation radius b, all tensor components can be brought 
above unity. If we start the transmutation from the boundary of the device (6 = 1), then 



A/" = 1 and we recover the results derived in 25 . 



A. 2 Invisible Sphere and Eaton lens 



Section 5. of 25 shows that transmutation is possible even when n(r) has a general form 
that diverges like as r — )■ 0. For example, such a general refractive index profile can 
be transmuted by the function R{r) = n{r')dr' , where UQ^b) is the normalisation 

constant. This transmutation will keep the angular eigenvalues of e*^ and constant and 
the radial eigenvalue will decrease slightly with increasing radius. For the Invisible Sphere 
the integral gives the expression 

= aTTm ^^^'^ /(O = Z , ■ (8) 

J^xsyo) n(r')dr' 

In the expression above the normalisation constant Mxs{b) was chosen such that R{b) = b 
and Mis{i) = 1- Therefore, Nxsip) = f{b)/b. In Fig. 3 we plotted Nxsip) in red as a 
function of b. 

Since Nxsip) increases as b decreases and n^, ne, and are all proportional to Afxs{b), it is 
easy to see that by a sufficiently small choice of b all tensor components can be raised above 
unity. Calculations show that b = 0.2489 is the largest and therefore optimal transmutation 
radius for which all the tensor eigenvalues of the Invisible Sphere are greater than or equal 
to one. In Fig. 4 we plotted the radial and angular eigenvalues for the partial transmutation 
with b = 0.2489. 
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Figure 3: The plot shows the normahsation constant for the Invisible Sphere (red) and 
the Eaton Lens (blue). Since Af increases as b decreases, by choosing a sufficiently small 
transmutation radius all entries of 5*- and fj,\ can be brought above one. 
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Figure 4: The plots show the values of the radial and angular elements of the mixed tensors 
e^j = fi\ of the partially transmuted Invisible Sphere, where b = 0.2489. Er changes between 
1 and 3.38, and ee and change between 1 and 11.38. The refractive index proffie remains 
unaffected outside the radius of the transmutation. 



For the Eaton Lens the integral gives 25 



i?(r) 



ATeib) 



f{r), where /(r) 



7r + 2 



arcsm^ / - + 



(9) 



In the expression above Nsib) ensures that R{b) = b and A/f(l) = 1. We plotted Meip) in 
Fig. 3 in blue. For the Eaton Lens the optimal transmutation radius is 6 = 0.3585. Fig. 5 
shows the radial and angular eigenvalues for the partial transmutation with b = 0.3585. 
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Figure 5: The plots show the values of the radial and angular elements of the mixed tensor 
e^j = fi\ of the partially transmuted Eaton Lens, where b=0.3585. Er changes between 1 
and 2.14, and ee and change between 1 and 4.58. 
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